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INTRODUCTION 

In the next decade ceramic composites will increasingly be used as high-temperature 
components in advanced gas turbines. Ceramic components will allow an increase in the operating 
temperature of the gas stream (measured by the turbine entry temperature), thereby resulting in 
greater fuel economy in aerospace and automotive turbine applications. Since these materials will 
be produced from abundant nonstrategic materials, computational structural mechanics methods 
are evolving to keep pace with this technology. As a result, establishing protocols for sound design 
methodology, which is the subject of this technical note, is the focus of much current analytical 
research. 

Our objective was to extend a reliability model proposed by Duffy and Manderscheid 
(1989) for orthotropic ceramic composites. Their approach used tensorial invariant theory to 
develop an integrity basis from which a subset of invariants was created to incorporate material 
anisotropy. Herein we propose a different subset of the original integrity basis and construct a 
more general scalar failure function. The assumption that the failure of a component is governed 
by its weakest link leads to a formulation that is similar in nature to the principle of independent 
action (PIA) model for monolithic ceramics. Note that this is a continuum approach to reliability 
analysis in that it excludes any detailed interaction between individual constituents. 
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EXTENSION OF THE RELIABILITY MODEL 

Consider a volume whose failure is assumed to be governed by its weakest link. Under 
this assumption, events leading to failure of a given link do not affect other links (see, for example, 
Batdorf and Heinisch (1978), Wetherhold (1983), and Cassenti (1984)); thus component reliability 
is given by the following expression: 


R = exp 



( 1 ) 


Here, Vfaj) is the failure function per unit volume at position Xj within the component, and V is the 
component volume. Note that the lower case Roman letter subscripts denote tensor indices with an 
implied range from 1 to 3. For orthotropic composites, the failure function must reflect the stress 
state and the appropriate material symmetry. This requires that 


aj.bj) 


( 2 ) 


where aj and b. are orthogonal unit vectors that represent the local principal material directions, 

and a- represents the Cauchy stress tensor. The sense of a^ and b^ is immaterial; thus their 

influence is taken through the products a- a- and b b-, that is 

^ J ^ J 


^oij,a i aj,b i bj) (3) 

Since the stress and local preferred directions may vary from point to point in the component, 
equation (3) implies that the stress field ^.(x^) and unit vector fields a-(x^) and b^x^) must be 
specified to define ip. 

Because ip is a scalar function, it must remain form invariant under arbitrary proper 
orthogonal transformations. Work by Reiner (1945), Rivlin and Smith (1969), Spencer (1971) and 
others demonstrated that by applying the Cayley-Hamilton theorem and the elementary properties 
of tensors, a finite set of invariants, known as an integrity basis, can be developed. Form 
invariance of ip is ensured if ip depends on invariants that constitute either the integrity basis or 
any subset thereof. We are not reporting the first practical application of this approach in 
formulating models that are functionally dependent on stress and material direction. Others who 
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have used this methodology are the following: Lance and Robinson (1971), who proposed a 
maximum shear stress plasticity theory for composites; Boehler and Sawczuk (1977), who proposed 
a plasticity theory for anisotropic cohesive materials; Arnold (1989), who developed a thermoelastic 
constitutive theory for transversely isotropic materials; and Robinson and Duffy (1990), who 
developed a viscoplastic constitutive theory for this same class of materials. Clearly, our direction 
here is not without precedent. 

Adapting the previously mentioned work to ip gives an integrity basis composed of 
twenty-eight tensor products. Following arguments similar to those of Spencer (1984), we find 
that several of these tensor products are equal and others are trivial identities such that the final 
integrity basis for ip contains only the invariants 


I = a- 
l ii 

(4) 

I = o-.o-. 
2 1J J1 

(5) 

\ = Vj^w 

(6) 

1. = Vf* 

(7) 

*s = a i a j°jk ff ki 
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( 10 ) 


From this group we can construct a subset of invariants that correspond to the stress components 
oriented to the principal material direction. This new set of invariants includes 


I = I 

(11) 
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For a uniformly stressed unit volume (or in the context of Weibull analysis, a single link) 

the first three invariants correspond to the magnitudes of the normal stress components. Note that 

I corresponds to the magnitude of the normal stress component in the direction defined by the 

cross product of vectors a. and b. , that is 

J * 


d. 
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6 ijk a j b k 


(17) 


where e^ is the permutation tensor. The last three invariants in the subset correspond to 
magnitudes of the shear stress components. (Refer to Fig. 1 for a pictorial representation of the 
invariants.) A fundamental difference between the theory proposed here and that of the original 
work by Duffy and Manderscheid (1989) lies in the number of invariants derived from the original 
integrity basis. The previous model employed only five invariants and suppressed the dependence 
of ip on invariants I and I . The current formulation uses six invariants and includes I , but it is 

2 3 2 

still independent of 1^. The additional parameter arises from a different and apparently more 
consistent set of assumptions. Taking 


i>= Hi ,i ,i ,i ,i ,i ) 
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( 18 ) 


ensures ip is form invariant. 

If we assume that compressive stresses associated with normals I 
contribute to failure, then 


I , and I do not 
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In addition, the shear stress contribution is sign-insensitive; thus 

<i,) = | i<| ( 22 ) 

(I 5 ) = | l5| ( 2 ^) 

<ia> = I ill ( 24 ) 


for all values of I , I , and I . At this point we assume that the stress components identified by 

4 5 6 

the aforementioned invariants act independently in producing failure (i.e. , a non-interactive 
theory). Following reasoning similar to that of Wetherhold (1983), we find that if) takes the form 


(I l) 1 

OH 

+ 

< I2> 1 

a 2 

+ 

' (I 3 > ' 

1 0i J 


_ P 2 . 


. P 3 . 


+ 

(i.)' 

a 4 

+ 

[<i 5 >l 

«5 

+ 

<i« > 1 


fit J 


1 

c n 

1 


/?• J 


(25) 


The a’s associated with each invariant correspond to the Weibull shape parameters, and the P’s 
correspond to Weibull scale parameters. Inserting equation (25) into the volume integration given 
by equation (1) along with equations (19) to (24) yields a reliability model for an orthotropic 
ceramic composite subject to a three-dimensional state of stress. 


CONCLUDING REMARKS 

In this technical note we have proposed a more general reliability model for ceramic 
composites with orthotropic material symmetry. This approach allows material orientation within 
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a component to vary along a family of curves such that the material is locally orthotropic. The 
model was constructed by using an invariant formulation. Such an approach indicates the 
maximum number and form of stress invariants necessary to define the failure function ip. A subset 
of the integrity basis for ip was constructed with invariants that correspond to the stress 
components in the principal material orientation. The result is a model similar in nature, yet 
different in form, to the PIA theory for monolithic ceramics. This model can be readily integrated 
with the finite element methods of structural analysis. Using the stress output from a finite 
element analysis allows component integrity to be computed by calculating element-by-element 
reliability. On the basis of the weakest link concept then, the component survivability is simply 
the product of individual element reliabilities. The model could be easily integrated into a 
computer algorithm such as SCARE (Structural Ceramics Analysis and Reliability Evaluation, 
Gyekenyesi (1986)). 
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Fig. 1 Physical interpretations of the invariants. 
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